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fH , 1. Introduction 

^ . Our ultimate goal is the construction of a model for interactions of two nu- 

Q I clei in the energy range between several tens of GeV up to several TeV per 

• nucleon in the centre-of-mass system. Such nuclear collisions are very com- 

plex, being composed of many components, and therefore some strategy is 
needed to construct a reliable model. The central point of our approach is 
the hypothesis, that the behavior of high energy interactions is universal 
5^ I (universality hypothesis). So, for example, the hadronization of partons in 

nuclear interactions follows the same rules as the one in electron-positron 
annihilation; the radiation of off-shell partons in nuclear collisions is based 
on the same principles as the one in deep inelastic scattering. We construct 
a model for nuclear interactions in a modular fashion. The individual mod- 
ules, based on the universality hypothesis, are identified as building blocks 
for more elementary interactions (like e~^e~ , lepton-proton), and can there- 
fore be studied in a much simpler context. With these building blocks under 
control, we can provide a quite reliable model for nucleus-nucleus scatter- 
ing, providing in particular very useful tests for the complicated numerical 
procedures using Monte Carlo techniques. 
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2. The Universality Hypothesis 

Generalizing proton-proton interactions, the structure of nucleus-nucleus 
scattering should be as follows: there are elementary inelastic interactions 
between individual nucleons, realized by partonic "half-ladders" , where the 
same nucleon may participate in several of these elementary interactions. 
Also elastic scatterings are possible, represented by parton ladders. Al- 
though such diagrams can be calculated in the framework of perturbative 
QCD, there are quite a few problems : important cut-offs have to be chosen, 
one has to choose the appropriate evolution variables, one may question the 
validity of the "leading logarithmic approximation" , the coupling of the par- 
ton ladder to the nucleon is not known, the hadronization procedure is not 
calculable from first principles and so on. So there are still many unknowns, 
and a more detailed study is needed. 

Our starting point is the universality-hypothesis, saying that the be- 
havior of high-energy interactions is universal. In this case all the details 
of nuclear interactions can be determined by studying simple systems in 
connection with using a modular structure for modeling nuclear scattering. 
One might think of proton-proton scattering representing a simple system, 
but this is already quite complicated considering the fact that we have 
in general already several elementary interactions . It would be desirable 
to study just one elementary interaction, which we refer to as "semihard 
Pomeron" , which will be done in the next section. 

3. The semihard Pomeron 

In order to investigate the semihard Pomeron, we turn to an even simpler 
system, namely lepton-nucleon scattering. A photon is exchanged between 
the lepton and a quark of the proton, where this quark represents the last 
one in a "cascade" of partons emitted from the nucleon. The squared dia- 
gram represents a parton ladder. In the leading logarithmic approximation 
(LLA) the virtualities of the partons are ordered such that the largest one is 
close to the photon [1, 2]. If we compare with proton-proton scattering, we 
have ordering from both sides with the largest virtuality in the middle, so in 
some sense the hadronic part of the lepton-proton diagram represents half 
of the elementary proton-proton diagram, and should therefore be stud- 
ied first. In fact such statements are to some extent commonly accepted, 
but not carried through rigorously in the sense that also for example the 
hadronization of these two processes is related. 

But first we investigate the so-called structure function F2, related to 
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the lepton-proton cross section via [3] 



da 2ti ( „ \ ^ , ^n, 

+ F2(x,Q2) (1) 
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and with the kinematic variables 

2pq 



Ft 

R=-Er-^, (2) 



Q' = -<l'; x = ^ (3) 



where q and p are the four-momenta of the photon the proton respectively. 
Fl is the longitudinal structure function. F2 represents the hadronic part of 
the diagram, and is, using eq. (1), measurable. In lowest order and consid- 
ering only leading logarithms of Q^, only two diagrams contribute , which 
turn out to be [2] 

j ^ ij ^0 ^ S S 

where P/ is the Altarelli-Parisi splitting function and Qo is some cut-off of 
the integration. The variable ^ is the momentum fraction of the quark 
with respect to the proton. Assuming some proton distribution /(^, Q^) at 
the "factorization scale" Qq) incoherent superposition provides [2] 

F2(x,Q2)=^e2xf (x,Q2) (5) 

3 

with 

PM') = f{xMi)+Y: £1 j^^%niM")^pi{^)- (6) 

Iterating this equation obviously represents a parton ladder with ordered 
virtualities . Strictly speaking, eq. (6) is still useless, because some of the 
functions P- (z) diverge for z ^ 1. However this is cured by considering 
virtual emissions, and this can be conveniently taken into account, by "reg- 
ularizing" the functions P/, which amounts to adding terms proportional 
S{l — z), which cures the divergence. For Monte Carlo applications it is more 
useful to proceed differently [4] : one distinguishes between "resolvable" and 
"unresolvable" emissions. Unresolvable emissions are the virtual ones and 
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emissions with very small momentum fraction (< e). Then one sums over 
unresolvable emissions which provides a factor 

A'(0§,0?)=exp{-i:/^:'f |"'<^«£^K)). (7) 

called Sudakov form factor. This can also be interpreted as probability of no 
resolvable emission between Qq and Qi- So the procedure amounts to only 
considering resolvable emissions, but to multiply each propagator with A*. 
Based on the above discussion, we define a so-called QCD evolution function 
Eqcb, representing the evolution of a parton cascade from scale Qq to Qf, 
as 

E^cBiQlQl^) = }Too^i^ciiQlQl^)^ (8) 

(n) 

where ^'qcd represents an ordered ladder with at most n ladder rungs. This 
is calculated iteratively based on 

sSoSCSo. e5.i) = s(i-x) Sij A-(o§, <5?) (9) 

+ E Q % r f i ^'-^^ ^^<i* 

The indices i, j, k represent parton flavors. The function -Eqcd is calculated 
initially for discrete values of the variables, and later used via interpola- 
tion. In this way we are really sure to use the same QCD evolution for 
any application, it is the same in deep inelastic scattering as in nuclear 
interactions. 

Next we have to determine the a;-distribution of the first parton of the 
ladder. We expect for the momentum share of this first (the fastest) parton 
of the ladder a distribution as l/x, which leads to mostly small values of 
X, which implies on the other hand a large mass M ~ l/x [5]. Such large 
mass objects are theoretically described in terms of Regge theory, the most 
prominent object at large masses being the Pomeron (P). So, as dicussed 
already in [6], the complete diagram is composed of the parton ladder, as 
discussed, and a soft Pomeron, given as 

^softP~X-°'^, (10) 

and the coupling between the soft Pomeron and the nucleon, which takes 
the form 

(11) 

We call -BsoftP also the soft evolution, to indicate that we consider this as 
simply a continuation of the QCD evolution, however, in a region where 
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perturbative techniques do not apply any more. We consider quarks and 
gluons to be emitted from the soft Pomeron, in case of quarks we have to 
spht the gluon momentum. So we have the following initial distribution, 

vip{^)=CTP^EU^, (12) 

with 

jpi _/ -E'softiP if ^ = 5 f-.o\ 

£;,oftP0P| iii = q,q ■ ^^^^ 

There is a second contribution, where a Reggeon is involved. The corre- 
sponding soft evolution is 

EsoftTR^-x"^^, (14) 

and for the Reggeon-nucleon coupling we take 

Cm ~ x-^^. (15) 
So we have the following initial distribution 

vi,{x)=Cu^EU^, (16) 

with 

^soft]R-|o ifi = g ■ 

The total initial distribution is obviously 

ip\x) = ipip{x) + ^iR{x), (18) 

which is by construction equal to P{x,Qq). The distribution at scale is 

i 

The structure function is then calculated as : 

F2(x,Q2) = ^e|x/^(x,Q2). (20) 

i 

For Q = Qo, the P-contribution is a function which peaks at very small 
values of x and then decreases monotonically towards zero for a:; = 1, the 
IR-contribution on the other hand has a maximum at large values of x and 
goes towards zero for small values of x. The precise form of / depends 
crucially on the exponent for the Pomeron-nucleon coupling, and we find a 
good agreement for /3f = i . 
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We are now in a position to write down the expression Ggemi ^oic a 
cut semihard Pomeron, representing an elementary inelastic interaction in 
pp scattering. We can divide the corresponding diagram into three parts . 
We have the process involving the highest parton virtuality in the middle, 
and the upper and lower part representing each an ordered parton ladder 
coupled to the nucleon. According to the universality hypothesis, the two 
latter parts are known from studying deep inelastic scattering, representing 
each the hadronic part of the DIS diagram 

p^J2f'^E^CD, (21) 



where 99* is the distribution at Q = Qq. Eq^^) represents the evolution 
between Qq and some Q^, is therefore a function of x and Q^. The 
complete diagram is therefore, for given impact parameter b and given 
energy squared s, given as 

Gsemi = E / dedCdQ^f\t,Q^)f{r,Q^)^^itrs,Q'), (22) 

ij ^ 

This may be written as 

Gsemi = E / dx-^dx- G'J^, {x+,x-), (23) 
IJ ■' 

with 

4Vmi(^+, = Ci{x+) Cj{x-) I du+du-dQ^ eUi ® 4cd(^+) 

ijkl 

EUj^ E'^CBi^-) -^{u+u-x+x-s, Q2). (24) 

The variables / and J may take the values P and JR. Integrating Ggemi 
over b, wc obtain the cross section to produce a pair of hard partons, which 
is measurable since the partons show up as hadron jets and may be re- 
constructed. This cross section is therefore called jet cross section, given 
as 

<^jet(s) = J dHGsemi{s,b). (25) 

So, in particular, checking the differential cross section da^et/dt against data 
provides an important consistency check. 

In addition to the semihard Pomeron, one has to consider the expres- 
sion representing the soft Pomeron [7]. The latter one, Ggoft, is the Fourier 
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transform of a Regge pole amplitude A ~ So an elementary inelastic 
intr action in an energy range of say 10 - 10^ GeV is therefore written as 

+ Gsoff (26) 



4. Hadron Production 

As discussed in the last chapter, there exist observables which can be cal- 
culated without detailed knowledge about hadron production. There exist, 
however, a huge variety of data concerning hadronic observables such as 
rapidity and pt spectra of different types of hadrons, multiplicity distribu- 
tions and so on. From a theoretical point of view, this requires some more 
details to be worked out. In addition, the numerical procedures are getting 
more complicated, and in fact, the most convenient method is provided by 
the Monte Carlo technique, which means in this context the following: for a 
given reaction, one writes the total cross section as a sum over contributions 
from different configurations X, 

<Ttot = E ^(^)' (27) 

such that o"(X) is non negative and can be interpreted as probability. The 
sum has to be replaced by an integral in case of a continuous configuration 
space. In general, such a configuration X represents a sum over a certain 
class of diagrams, because individual diagrams may be negative. So our 
work is, besides providing appropriate expressions for the total cross sec- 
tions, also to provide algorithms to generate configurations X according to 
the corresponding distributions. 

The first step amounts to generating parton configurations before wor- 
rying about hadronisation. 

Let us start again with the case of lepton-nucleon scattering . We first 
have to generate the kinematical variables x and according to the lepton- 
proton cross section, 

- ' " - 2y + F2{x, Q2). (28) 



Then a Pomeron or Reggeon type of coupling is chosen, with probabilities 
/F2 and F2^/F2 respectively. We have 

F^/'^ = xY,e]vl>/^^^CB- (29) 

ij 
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Eqcd can be written as 

^§cd(Qo, z) = 6{1 - z) 6ij A\Ql + &^^^{Ql q\ z), (30) 

where the first term represents no parton emission and the second term at 
least one emission. The probabihty of no emission is therefore 

--^xY.e]^]p/^{x)Ai{QlQ% (31) 

^2 3 

and a configuration with no emission is generated according to this proba- 
bihty. In case of at least one emission, one generates the flavour i and the 
fraction xq of the first parton of the QCD cascade according to 

E V V^Wh(^o) ^gcD(Qo, Q\ f )• (32) 

j XQ XQ 

We arc left with the problem of generating the cascade of partons starting 
from a parton with fraction xq, flavour i, at scale Q^, up to the photon 
vertex at scale Q"^. We have 

E%CT^{QlQ^^)= (33) 

and therefore k, z\, and Q\ of the next emission are generated according to 
^,2Ai^lg^,,^^,^.^_^(gj^j,^,_ (34) 

and so on. This completed the description of the algorithm to generate 
parton configurations, based on exactly the same formulas used to calculate 
F2 as discussed in the previous section. 

The next step consists of generating with certain probabilities hadron 
configurations, starting from a given parton configuration. We cannot cal- 
culate those probabilities within QCD, so we simply provide a recipe, the 
so-called string model. The first step consists of mapping a partonic config- 
uration into a string configuration. For this purpose, we use the colour rep- 
resentation of the parton configuration: a quark is represented by a colour 
line, a gluon by a colour-anticolour pair . One then follows the colour flow 
starting from a quark via gluons, as intermediate steps, till one finds an 
antiquark. The corresponding sequences 



q- gi- 92- ■ ■ - 9n- q 



(35) 
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are identified with kinky strings, where the gluons represent the kinks. 
Such a string decays into hadron configurations with the corresponding 
probabihties given in the framework of the theory of classical relativistic 
strings. 

The above discussion of how to generate parton and hadron configura- 
tions is not yet complete : the emitted partons are in general off-shell and 
can therefore radiate further partons. This so called timclikc radiation is 
taken into account using standard techniques. The mapping of parton to 
hadron configurations still works the same way as discussed above. 

Let us now discuss the generation of parton configurations for an el- 
ementary proton-proton interaction represented by a semihard Pomeron. 
Based on 

Gsemi = J2 dx^dx~ (36) 
IJ 

one generates the light come momentum fractions x"*", x~ of the "Pomeron 
ends" according to G. This quantity may be written as 

Gli^{x+,x-) = Cj{x+)Cj{x-) (37) 

■ dz+dz- EU,{z+) El,,j{z-) 4^z-x^x-s), 



^3 



with 

a]4(s) = J2 f dw+dw-dQ^ (38) 

kl '' 

flfykl 

The integrand serves as probability distribution to generate z'^ and z~ . 
The two quantities z"*" and z~ are the momentum fractions of the "lad- 
der ends" with respect to the "Pomeron ends". Knowing the ladder mass 
s = z~^z~x'^x~s, we have to generate the ladder rungs. Generalizing the 
definition of Ujct, we define 



kl 



Ei^CDiQl Q\ ^'^) E'^cb{QIq\ W-) -^{w+w-s, Q2) 



Ql^) = JlJ dw+dw-dQ^ (39) 

'kl 
dQ^ 

<id(Q?'Q2,^) = Y.ldw-dQ' (40) 

dQ2 



and 

i2 s\ 



f]~-ki 

eScMQ', W-) A\Ql Q2) ^^{w-s, Q'). 
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representing ladders with ordering of virtualies on both sides (<Tjet) or on 
one side only (cTord)- We calculate and tabulate crjet and (Tord initially so 
that wc can use them via interpolation to generate partons. The generation 
of partons is done in an iterative fashion based on the following ladder 
equation: 

^UQlQl')=<rniQlQls) (41) 

k 

Our treatment of generating parton configurations for an elementary 
pp interaction is absolutly compatible with deep inelastic scattering, it is 
based on the same building blocks, in particular on the evolution functions. 

Our procedure can now be used to treat the case of many semihard 
Pomcrons, for p — p as well as nuclear scattering. A detailed discussion will 
be given in a future publication. 

This work has been funded in part by the IN2P3/CNRS (PICS 580) and 
the Russian Foundation of Fundamental Research (RFFI-98-02-22024). We 
thank the Institute for Nuclear Theory at the University of Washington for 
its hospitality and for partial support during the completion of this work. 

References 

1. E. Reya, Physics Reports 69 (1981) 195 

2. G. Altarelli, Physics Reports 81 (1982) 1 

3. R.K. Ellis, W.J. Stirling, and B.R. Webber, QCD and Collider Physics, Cam- 
bridge University Press, 1996 

4. G. Marchesini and B.R. Webber, Nucl. Phys. B238 (1984) 1 

5. A. Donnachie, P.V. Landshoff, Phys. Lett. B332 (1994) 433 

6. K. Werner, H.J. Drescher, E. Purler, M. Hladik, S. Ostapchenko, in proc. of 
the "3rd International Conference on Physics and Astrophysics of Quark-Gluon 
Plasma", Jaipur, India, March 17-21, 1997 

7. K. Werner, Physics Reports 232 (1993) 87 



